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Preface 


This note serves as an elementary introduction to the induced velocity created by a field of 
vortices that reside in the wake of a rotor blade. The approach is to build a bridge between familiar 
fixed-wing theory over to a rotor blade in forward flight. This bridge is built in four parts plus 
concluding remarks. 

By way of background, rotorcraft technologists, after a nearly seven decade effort, have 
finally begun to provide accurate computer based prediction of a helicopter rotor blade’s airloads 
and dynamic response. The decades-long effort was capped during the last four years by tying 
airload prediction (with advanced CFD methods) to completely coupled structural dynamic response 
(calculated with very advanced modal methods). Typical results, such as Refs. 1 and 2, confirm the 
progress of what this relatively small band of engineers have achieved after some 70 years of 
dedicated work. 

In September 1969, 1 presented a paper at the V/STOL Technology and Planning Conference 
sponsored by the Air Force Flight Dynamics Laboratory (Ref. 3). 1 This paper included a figure 
showing the progress in removing assumptions from the original rotor performance theory 
developed by Juan de la Cierva in the late 1920s. The figure is reproduced in this note as Figure 1. 
You can see from Fig. l’s dashed line that I was, in 1969, confident that by 1975 we would be done. 
Of course, a revision to that 1969 view reflecting history is quite in order. So, Fig.l’s solid line now 
shows the more accurate progress in hindsight. 

One of the toughest assumptions to remove was that of uniform downwash. Professor Rene 
Miller of MIT showed the way in 1962 with Ref. 4. Work that followed is discussed rather 
completely in Chapter 13 of Wayne Johnson’s Helicopter Theory. The interaction between the rotor 
blades that create the induced velocity field and the effect the induced velocity feeds back on each 
blade was a daunting complexity that needed a computer before even rudimental solutions became 
tractable. 

There is some interesting knowledge to be gained, however, from disconnecting the blade 
loading from the induced field. That is the purpose of this note: to calculate the induced velocity at a 
blade where the blade’s bound circulation is given in terms of radius and azimuth. Once specified, 
the blade’s bound circulation and lift distribution remain unchanged despite the resulting induced 
velocity field. This is, of course, a comparatively simple problem compared to the real problem. 
Furthermore, in this note only a one-bladed rotor is considered, and to make the problem even 
simpler, only a prescribed, rigid wake is considered. 

An interesting advantage of these simplifications is that a bridge from the fixed wing to the 
rotary wing is quite easy to construct. To make the bridge secure, I have tried to leave nothing to the 
imagination concerning the sign conventions, the mathematical notations, and the steps under 
discussion. 

F. D. Harris, Oct. 16, 2005 


1 This conference was held in Las Vegas, Nevada, a city — I was led to believe — that was the V/STOL technical capital 
of the world. 
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I. Introduction to the Fixed Wing 


The Classical Fixed-Wing Problem 


The classical fixed-wing induced velocity problem begins, as explained in most aerodynamic 
textbooks, by assuming a bound circulation that is elliptical. The wing is placed in straight and level 
flight. An array of horseshoe vortices is envisioned. The resulting induced velocity at the wing is 
found to be constant from the port wingtip to the starboard wingtip. In summary, 


(0.1) if 


= r 



then 


2b 


where T 0 is the maximum circulation in square feet per second, b is the wingspan in feet, and xw is 
the span station being xw = - b/2 at the port wingtip and xw = +b/2 at the starboard wingtip. The 
induced velocity, v xw , is in feet per second and is constant from tip to tip. 


In this note, think of the wing rotating, not flying straight. Therefore it is quite helpful to 
approach the fundamental geometry used in the classical derivation from a different point of view. 
Suppose a fixed wing is east of a pylon and flying north past the pylon. Ignore — for the moment — an 
anticipated 180-degree U-turn around the pylon.' (The U-tum results will be discussed after 
reconstructing the classical problem.) The situation is illustrated with Fig. 2. 


In Fig. 2, the right-hand axis system gives a positive Z axis coming up out of the paper. The 
“pylon” is located atX = Y = Z = 0. The wing centerline is placed a distance D from the pylon in the 
plus X direction. The wing coordinates are measured from the wing centerline. The station on the 
wing where induced velocity is sought is denoted by xw. A vortex trails downstream from the wing 
at wing station xv. A small element of the vortex, dS, is shown located at a distance L from wing 
station xw. The reference angle 9 is used to locate the vortex element in relation to the wing. When 
9 is equal to minus tc/ 2, the vortex element is at -oo. When 9 equals zero, the vortex element is 
located at the wing. The calculation of induced velocity by the Biot-Savart law, as derived from 
vector notation, is simply 


( 0 . 2 ) 


d(dv xw ) 


Yv 

i 

r 

"q. 

jy) 

1 

r 

"q. 

yi 

i 

de -_L 


f dr 2 

xw 

dxv 

1 

J 


LjdSj - LjdSj 

4n 

L 3 

4n 


^ dxw j 

XW— XV 

L 3 


This fundamental equation is deceptively simple because to calculate the induced velocity, v x 
any station along the wing, xw, only a double integral has to be performed. That is 


at 


(0.3) 


1 

/»+b/2 

0 

~ 

( dr Y 

xw 

dxv 

1 


LjdSj - LjdSj 

47T 

J 

-b/2 J 

-Till 


^ dxw j 

xw=xv 

L 3 


dO 


2 Imagine a rotor blade at the traditional downwind, zero-azimuth station, which will, after the U-turn, be at the 180- 
degree upwind azimuth station. 
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The double integration required by Eq. (0.3) is hampered (to put it mildly) whenever — — 

dxw 

equals either plus or minus infinity. The double integration can be an even bigger problem whenever 
L = 0. If the double integration is performed numerically, situations where numbers approach plus or 
minus 10 10 can become quite frustrating. Modem computers using so-called double precision and 
advanced numerical integration schemes have helped to lower this frustration. But the fundamental 
basis of numerical integration is the Taylor series. The hope is that — without a breach in engineering 
accuracy — d9 and dxw can be replaced by A9 and Axw when the integral operators are replaced by 
summation operators. There are, of course, a number of ingenious coordinate transformations that 
can completely remove an apparent integrating roadblock. 

Here is an example of a coordinate transformation that helps lower frustrations with the 
fixed-wing problem. Consider the situation when a wing’s lifting line is loaded with an elliptical 
bound circulation defined as 


(0.4) 


r 


= r 



where xw = - b/2 at the port wingtip and xw = + b/2 at the starboard wingtip. The derivative that the 
Biot-Savart law requires is 


(0.5) 


Yv = 


-( ^ xw dxw ) 
V dxw J 





K 





dxv 


where xv = -b/2 at the port wingtip and xv = +b/2 at the starboard wingtip. You can immediately see 
that the required derivative is plus infinity at the starboard wingtip and minus infinity at the port 
wingtip. Now look what happens with the coordinate transformation of 

b b 

(0.6) xv = - — cosP and dxv = — sin P dp 

where P = 0 is the port wingtip and P = n is the starboard wingtip. Then the bound circulation 
becomes 


(0-7) r 3 =r l- 


r 

f 2) 

2 

f b 

J l- 



cosP 


IbJ 


2 ) 


cosP =T o ^/ 1-cos 2 P =T o sinp 


and the required Biot-Savart derivative becomes 


( 0 . 8 ) Yv =- 


dT 


dxw 


-dxv 


\ 

J XW— XV 


-y/l-COs'’ p 


( b ) 

(— cos P) —sinp dp 

\2 J 


= -T 0 cos P dp 
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Clearly, singularities caused by the derivative — — at the wingtips have been removed. 

dxw 

The preceding coordinate transfonnation, when placed in the Biot-Savart law for this fixed- 
wing problem, gives 


r K r° 


(0.9) 


v = 

xw 


4n 


[-r o cos p dp] 


0 J 


-n/2 


LjdSj-LjdS; 

L 3 


d9 


The next step is to construct the second portion of the integrand, which is 


LjdSj -LjdSj 

L 3 


de 


The geometric dimensions, L, L,, Lj, dS, and dSj , are components of vectors. When looking at 
Fig. 2, one has a choice of reference systems. Since this fixed wing is flying straight and level past 
the pylon, there is no advantage to including the distance D or using the reference angle 9 in the 
wake geometry. 3 Therefore, continue using just the reference angle 8. That means 


MSj-LjdS; 

L 3 


d9 is replaced by 


LjdSj -LjdS; 

[} 


dS 


and 8 is positive counterclockwise. The vortex segment farthest downstream is located at 8 = 9 and a 
vortex segment right at the wing is located by 8 = +n/2. 


Now, from Fig. 2, the distance L is written as 

L = Lj +Lj +L k = (Xw -Xv)i + (Yw- Yv)j + (Zw-Zv)k 
with Xw = xw Yw = 9 


(9.19) 


Xv = xv 


Yv = 


xv -xw 
tan 8 


and L 3 = [(Xw-Xv) 2 +(Yw-Yv) 3 +(Zw-Zv) 2 ] 
Therefore, with the usual “substitute and simplify” phrase, 
(xv-xw) 


3/2 


(9.11) L = 


sin 2 8 


3/2 


Zw = 9 
Zv = 9 


In a similar manner, the vortex is described by the vector 
(9.12) S = S 1 +S J +S k =(Xv)i + (Yv)j + (Zv)k 


3 The next part of this note addresses the U-turning wing. That problem requires, of course, some reference system that 
includes D and uses the reference angle 0. 
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and therefore 


dS,=^dS = ^d8 = 0 

1 as as 


(0.13) dSj =— ^-dS = 

as 


aYv (xv-xw) 

dS = - - 

aS sin - 5 


, c as k azv 

dS. = — -dS = d8 = 0 


as 


as 


Here the assumption about the trailed vortex geometry is that it extends straight aft of the wing 
without descending or climbing, remaining perpendicular to the wing lifting line. Again, by 
substitution and simplification, 


(0.14) 


LjdSj -LjdSj 

L 3 


dS = 


-1 


xv-xw 


-sinSdS 


The double integral created by the Biot-Savart law now appears, after substituting Eq. (0.14) 
into Eq. (0.9), as 


(0.15) 


4n 


r* 

-» nil 

[-r o cospdp] 

-1 . s " 
sinS 

Jo*' 

0 

xv-xw 


dS 


A minor problem is immediately observed in Eq. (0.15). Fonnal integration really requires that xv 
and xw be related, in some fashion, to P and/or 8. This minor problem is repaired using Eq. (0.6) 

where xv = -^-cosP . And, rather than have the vortex wing station keyed by the angle P while the 

wing station is keyed to xv, you can make the substitution thatxw = -^-cosa . These substitutions 
result in 


(• 11 /- ti /2 


(0.16) 


v„ = ■ 


2tm 


[-r o cos p dp] 


0 J 


-1 

cos a- cos P 


-sin 5 


dS 


A much, much bigger problem with Eq. (0.15) and Eq. (0.16) is that any integration, whether fonnal 
or numerical, must face the possibility that (xv - xw) might be zero. In Eq. (0.16) terms, real trouble 
will occur if P = a. 


The integration with respect to wake age (i.e., with respect to S ) results in 


(0.17) 


F o 

f* 71 

cosp 

dB - 

/* 71 

cosp 

2b7T ^ 

0 

cos a- cos P 

P 21)71 ^ 

0 

1 

O 

0 

TO 

1 

O 

0 

c n 

P 

1 


Eo. 

2b 


Note immediately that the above derivation gives a negative value for induced velocity created by a 
wing carrying positive lift! The reason for this outcome is the right-hand rule axis system of Fig. 2 
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where the Z axis is positive up. Therefore, the induced velocity is directed in the negative Z 
direction. Of course, reason prevails in all textbooks on the subject, so the negative induced velocity 
is simply called downwash and the sign is changed to positive. This is equivalent to letting 
Yv = + r o cos p dp. 


Incidentally, Glauert proved the integration with respect to P equals n long ago. Also, Alan 
Pope in Appendix 3 of his book Basic Wing and Airfoil Theory shows that 




v„ = ■ 


(0.18) 


2b7t 

2lm 


cosp 


cos p- cos a 


r 

dpn — — cosa 
o 21m 


r 

S' 71 ~ 

dp= » 


2b7r 

J J 

0 - 


1 + - 


cosa 


cos p- cos a 


dp 


1 


cos p- cos a 


dp = [ ttI + -^- 2 - cos a [ 0] 

2b7r L J 2b7r L J 


It comes in very handy at times to know that 


(0.19) 



1 

^ 0 

cos p- cos a 


dp = 0 


and that 


pn 

cosnP 

An 7rsinnP 



_ . Q 

^ 0 

cosp -cos a 

sinp 


But now suppose numerical integration (rather than formal, closed-form integration) of Eq. 
(0.16) is the solution approach. I took an approach with the Mathsoft Inc. product called Mathcad 
Plus 6.0 that goes like this: 


<N 

II 

o 

u 

b = 

1 

8 = 0,1.... 

...179 


a s =(s + 

oo s 
O 


dvdpd8 n 

_ r o 

cosp m 

,m 2b7T [ 

cosP m -cosa s 

n = 0,1... 

90 

o„ = n 

n 180 

m = 0,1.. 

....180 

o _ 

~ 180 m 


dvdP - 71 y dvd P d5 n. m +dvdpd8 n+L 
Pm 180^ 2 

v 7T y dvdp m +dvdp ni l 
s 180^ 2 


b 

XW C = COSOL 

s 2 


This numerical solution, “programmed” in Mathcad, gave +0.9999746151 accuracy nearly within a 
blink of an eye, using a Dell Optiplex high-end computer. The answer from Eq. (0.17) is exactly 
unity (when you ignore the minus sign). 
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As a reminder, the fixed wing represented by a lifting line and an elliptical bound circulation 
produces a lift calculated from 


( 0 . 20 ) 


L = pvr 0 


^+b/2 

i 

-b/2 

Till 


1- 


^2xw^ 2 


= pvr r 2 a/ 1 -sin 2 0 

0 J - jc/2 


b j 

f 


dxw 


— cosOdO 
2 


b f+tt/2 ^ 

= pVT — f cos‘0d0 

0 2 J-n/2 


= pvr 0 

0 2 1 2 


= - P bvr 0 


so 


4L 

Ttpb V 


The wing incurs an induced drag to carry this lift, which is 


v ( n V r 

(0.21) Induced drag = L^^= — pbVT — — 

V U Jy2bV 


\ n jJ 

y~8 pr ° _ (p/2)V 2 7ib 2 


This drag can be used to calculate a horsepower required by multiplying both sides of the equation 
by velocity, V, and dividing by 550. Thus, 


( 0 . 22 ) 


Induced Horsepower = HP 


V 

L 2 1 

. Vq 

f L ) 

550 

i 

To 

rT 

CL 

1 

£ 

o 

IT) 

IT) 

1 

Ub 2 J 


w itli q = (p/2)V 2 


A numerical example calculated here will be a useful result to which one can later compare 
rotating-wing calculations. Suppose both a wing and a rotor have equal span, say, 44 feet. In rotor 
notation, the rotor radius is 22 feet. Assume, for example, that each lifting device carries a lift of 
2,712 pounds, each is flying at sea level (p = 0.002378 slug/ft 2 ), and each is flying at V = 301.8 
ft/sec or 178 knots. Then, by Eq. (0.22), the induced horsepower required by the wing is 6.13 
horsepower. 


Keep in mind that the elliptical bound circulation distribution used in the above discussion is 
the first term of the more general distribution used in fixed-wing analyses. Recall that the general 
distribution is seen in the form 

00 

(0.23) T = 2bV^A n sinnp 

1 

Everything read in this note — including the rotor analyses — could be extended by this fixed-wing 
logic o/Eq. (0.23). 


And there you have the fixed wing flying north past the pylon located off its port wingtip. 
Now consider the situation after the fixed wing does a 180-degree U-turn. 
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The Fixed Wing after Completing a 180-Degree U-Turn 


The next step in bridging the gap between a fixed wing and a rotary wing was actually taken 
by H. Glauert in 1923, although it is doubtful he had a bridge in mind when he published R & M 
866. After all, in 1923, Juan de la Cierva’s earliest autogyro experiments were just bearing fruit in 
Spain and a practical helicopter was still 15 years away. The title of Glauert’s 1923 report is 
“Calculation of the Rotary Derivatives Due to Yawing for a Monoplane Wing.” He was dealing with 
the wing rolling and yawing moments created by an induced velocity field that trailed the wing in 
one half of a circle. The problem is, in fact, quite akin to a hovering one-bladed rotor where only a 
small part of the wake is taken into account. Glauert obtained a very simple closed-fonn solution by 
assuming the wing made the U-turn with a large turning radius, D, relative to the wingspan, b, (i.e., 
D»b). His quite useable engineering result for the induced velocity over the wingspan was (in the 
notations of my note here) 


(0.24) 


lofl + I^ 

2b( 2D, 


The purpose of this portion of the note is to explore this problem when the turning radius is 
considerably smaller than what Glauert assumed. The objective is to think of the “pylon” as a rotor 
hub and place the port wingtip a small distance from the hub, which gives the appearance of a “root 
cutout” in rotorcraft tenninology. Glauert’s approach will be presented first and then a numerical 
integrating approach will be shown. 


The geometry of the problem is illustrated with Fig. 3. In contrast to Fig. 2, the wing has now 
advanced 180 degrees. Immediately note that the X axis is positive to the left, but left the Y axis in 
its commonly found, ordinate position. By the righthand rule then, the positive Z axis now points 
down, which is into the paper. The wing is doing a U-turn of distance D, which is measured from the 
“pylon” to the wing mid-span point. The wing is represented by a lifting-line vortex having an 
elliptical bound circulation, just as with the classical fixed-wing problem. The Biot-Savart law is 
again invoked, so 


(0.25) 


d(dv xw ) 


v L.dSj -L.dSj 

i v 1 J J 1 

4n L 3 


de 


and the dimensions are expressed as vectors. However, the basic vortex geometry has changed in 
going from Fig. 2 to Fig. 3. Now 


L = L; + Lj + L k =(Xw-Xv)i + (Yw-Yv)j + (Zw-Zv)k 

Xw = D + xw Yw = 0 Zw = 0 

Xv = (D + xv)cos0 Yv = -(D + xv)sin0 Zv = 0 

L 3 = [(Xw - Xv) 2 + (Yw - Yv) 3 + (Zw - Zv) 2 ] 3/2 

The reference angle 0 is taken positive counterclockwise in Fig. 3. In a similar manner, the vortex is 
described by 


where 

(0.26) 

and 
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(0.27) S = S,+S J +S k =(Xv)i + (Yv)j + (Zv)k 

but with changed vortex geometry reflecting the V 2 -circle wake, the vortex element geometry is 

5S ; 5Xv 5(D + xv)cos9 , 

dS, = — L d9 = d9 = — d0 = -(D + xv)sin0 

1 dQ dQ dQ v ’ 

dYv 5f-(D + xv)lsin0 

(0.28) dS ; =^d9 = d0 = — d0 = -(D + xv)cos0 

J dQ dQ dQ v ’ 

d s. =-^de=— de=o 


dQ 


dQ 


Here the assumption about the trailed vortex geometry is that it extends in a circular arc aft of the 
wing for 180 degrees without descending or climbing. Furthermore, any given vortex has constant 
circulation, y v , from when it leaves the wing all the way back to when the turn began (i.e., 9 = 0 
back to 9 = - 7t). 


The substitution of this U-turn geometry into the Biot-Savart law gives, with simplification, 


(0.29) d(dv„) = )y 

471 


(D + xv) (D + xw ) cos 9 - (D + xv)" 


(D + xv) 2 +(D + xw)" -2(D + xv)(D + xw)cos0 


— i3/2 


d0 


The wake-age integral problem is immediately seen as requiring elliptical integrals — even if 
the wing’s bound vortex circulation varies with wake age. For this example, assume as Glauert did, 
that y v does not vary with 9. Now, at the risk of boring you, the transfonnation to complete elliptical 
integrals (i.e., E and K) begins by letting 9 = k - 2cJ) and d9 = - 2 dc(). You also need to recall that 

cos 9 = cos (71 — 2(j) j = -cos(2())) = 2sin 2 cj)-l 
Then, a couple of substitutions and rearrangements immediately shows that 


d(dv xw ) = 


(0.30) 


y v 2(D + xv)(2D + xw + xv) 

1 

4tt (2D + xw + xv) 3/2 

|^1 — k 2 sin 2 

*] J 

y v 4(D + xv)(D + xw) 


sin 2 (j) 

=■ dc() 

4tt (2D + xw + xv) 3/2 

Too — 13/2 

1-k sin (j) 


d4> — 


where 


(0.31) k : 


4(D + xw)(D + xv) 
(2D + xv + xw)" 


This substitution changes the wake-age integrating limits from 9 = -n to 9 = 0 over to (j) = 0 
to (j) = +7t/2. The two budding integrals can be found, for example, in the translated Russian 
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handbook by I. S. Gradshteyn and I. M. Ryzhik titled Tables of Integrals, Series, and Products, 
edited by Alan Jeffrey. Thus 


(•n/2 

0 

/»7t/2 


(0.32) 


|^1 — k 2 sin 2 (j)j 

sin 2 (j) 
1-k 2 sin 2 (j)j 


3/2 # = 
3/2# = 


1 


1-k 2 

/» tt/2 


f Till 


Too -|3/2 

1-k" sin" (j) 


dcj) — 


COS 2 (j) 


to n3/2 

1-k" sin" (j) 


d(j) 


1-k 2 


E- 


-^-K ^E 


where the complete elliptical integrals, E and K, are computed as 


run 

1 d<t> 

-k 2 sin 2 (j) 

and their values depend on the modulus, k 2 , which for this U-turning wing problem, is given by Eq. 
(0.3 1). Note that when xv = xw, k 2 = 1 .0, E = 1 .0, and K = +oo. 


(0.33) E= f Vl-k 2 sin 

J 0 


K = 


In this way 


Yv 




(0.34) dv xw = f- 
4n 


(D + xv) (D + xw) cos 9 - (D + xv)‘ 


(D + xv)~ +(D + xw) 2 -2(D + xv)(D + xw)cos9 


-i3/2 


dO 


is converted into 


dv v , , = 


(0.35) 


y v 2(D + xv)(2D + xw + xv) 
4n (2D + xw + xv) 3/2 

y v 4(D + xv)(D + xw) 


(• Till 


3/2 


47T (2D + xw + xv) 
which, upon simplification, reduces to 


rn/l 


T 

< 


sin 2 (j) 



o o “|3/2 


r 

1-k" sin" (j) 

J 0 

IL 

-1 


|^1 — k 2 sin 2 (j)J 
dcj) 


d4> — 


(0.36) dv„ = f 
4 k 


K 


1 — 

xv - xw 2D + xv + xw 


The wake integration being complete, the problem is reduced to the spanwise collection of all 
vortices trailed from the wing. 


The insightful step Glauert next took was to make use of the approximations for E and K 
when k is closer to 1.0 than to 0 rather than calculate them using Eq. (0.33). These approximations 
are 
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E*l + 


(0.37) 


ifhrf^UVk 2 ) 

4 v U-k 2 J 2) y ’ 


v i, r i 6 i i , f i 6 

2 1 1 -k J 8 ll-k 2 


i (i-k 2 ) 


The fact that the approximation for E begins with one (1) is enormously important, which is seen 
when these approximations are substituted into Eq. (0.36). Of course, a lengthy integration problem 
results, which is 


(0.38) 


v = 

xw 


4 71 V xv - XW 


( — L ~) 1 ln ( -^rV-lfi-k 2 ) 

47rl xv -xw ) 4 1 1 1 — k J 2 f ’ 


Y v | 1 ] 1J 16 

4n v 2D + xv + xw J 2 vl-k 2 


Yv f 1 ] 1 J 16 

4ti: 1 2D + xv + xw j 8 1 U-k 2 


i (i-k 2 ) 


But now look very closely at the first integral to be obtained. Recognize that 


v = 

xw 


Yv [ 1 


4n v xv - xw 


is nothing more than the classical fixed-wing problem presented in Part I of this note. 

This is a key result that Glauert provided in his 1923 report because it says the so-called 
“near wake” of a U-turning fixed wing (i.e., think a rotor blade’s near wake) is no tougher problem 
than the straight flying wing to which is added the influence of a curved “far wake.” For the 
elliptical bound circulation used as the example in this note, the spanwise integration becomes 


(0.39) 


47tlxv-xw) 4l vl-k 


I L JE )_l (l _ k ») 

A 1 _U 2 0 V > 


4n v 2D + xv + xw J [_ 2 v 1 ~ k , 

i ifjTf 

4n l 2D + xv + xw ) 8 1 l 1-k 


i (i-k’) 
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The additional three integrals can be grouped into one integral. But first, the vortex 
circulation strength, assuming an elliptical bound circulation for the wing lifting line, Eq. (0.8), is 
substituted into Eq. (0.39). Next the elliptic integral modulus, k~ , from Eq. (0.31), is substituted into 
Eq. (0.39). (I would have included the results of these substitutions, but the resulting expression is 
way too long.) Then a selection of wingspan stations, xw, is made (say 50) and, in short order, 
Mathcad calculated the induced velocity distributions at any turning distance, D. 

Despite the appearance of possible singularities in Eq. (0.39), Mathcad actually experienced 
no numerical integration problems. I did not let Mathcad try to simplify the integrand. The 
temptation is to fiddle with 1-k , which Mathcad or I fouled up. A very unproductive effort given 
Mathcad’ s speed on my Dell computer. 

Glauert, in R & M 866, assumed that the turning radius was considerably greater than the 
wingspan, which allowed simple integration of Eq. (0.39), the result being Eq. (0.24). But consider 
the results as the wing makes tighter and tighter U-tums. A non-dimensional measure of the 
semicircle’s tightness is wingspan, b, divided by the distance from the “pylon,” D. The tightest turn 
would be when the port wingtip is touching the “pylon,” in which case D = b/2 or b/D = 2.0. This 
corresponds to a rotor blade with zero root cutout. A turning ratio of b/D = 1 corresponds to a 0.33 
root cutout. The extreme in the other direction would be, of course, not turning at all and so D = oo 
and b/D = 0. 

For simple illustration purposes, let the wingspan be unity (i.e., b = 1.0 foot) and let the 
maximum elliptical bound circulation, which occurs at the wing mid-span, be two (i.e., T 0 = 2.0 
square feet per second). Thus, for the following examination, T 0 /2b =1.0 foot per second. Two 
results are immediately known from the preceding discussion: 


Turn parameter, 
b/D 

Spanwise induced 
velocity, v x 

Source 

Not Turning 

v x = T 0 /2b= 1.0 

Eq. (0.1) 

Wide Turn 

v x = l+xw/2D 

Eq. (0.24) 


The results for several tighter and tighter turns are shown with Figures 4 through 7. For a 
turn distance D that is 10 times the wingspan, Glauert’s approximation is very useful as Fig. 4 
shows. When the turn distance is equal to the wingspan (equivalent to a rotor blade root cutout of 
1/3), the distortion in induced velocity across the wingspan is significant, as shown with Fig. 5. This 
distortion grows more pronounced as the turn distance shrinks, which is illustrated by Figs. 6 and 7. 
These two figures correspond to a rotor blade root cutout of 0. 1667 and 0.0476, respectively. 

Now consider the numerical double integration involved in this U-turning wing problem. 
Figure 6 offers a virtually exact example to which numerical integration can be compared. The 
integration problem at hand is this: 
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(0-40) v xw = 


/»+b/2 

f0 r 


>1 

| 


4n 

-b/2 J 

-n 


(D + xv)(D + xw)eos0-(D + xv)*' 


(D + xv) +(D + xw) -2(D + xv)(D + xw)cos9 


3/2 


dOdxv 


Philosophically, Glauert found the fixed-wing induced velocity equation is buried obscurely within 
this U-tuming wing problem. Therefore, the Mathcad-provided integration scheme will encounter all 
of the fixed-wing problems at the wingtips. (Mathcad returns the notice “will not converge,” which 
makes for frustration.) A second thing you kn ow is that the integrand’s denominator can get very, 
very small or even be zero when 9 = 0. Both these probable problems suggest using the fixed-wing 
coordinate transformations, which to repeat, are 


(0.41) xv = -— cosp 


dxv = — sinB 
2 


xw = — cos a 
2 


and assume an elliptical bound circulation along the wing’s lifting line so that 
(0.42) y v =-r o cospdp 

With the above thoughts in mind, a “brute force” numerical integration goes like this: 


r o =2 b = 1 D =0.7 

Dimension integration. 

Number of radial stations at which vortices leave wing, M = 90 
Number of azimuthal stations between 0 and 180 degrees, N = 18,000 
Range of radial stations where induced velocity is calculated, s = 0, 1 M - 1 


Then proceed with these calculations 


1 

a = — I s + — 


p -- | s) m 


M 

n = 0,1 N 

m = 0,1 M 


r 

■ — ( 
4tc 


( D 

-| cosP -) 

(D-^cosa s jcos0„-^D-^cos(3 m 

2 

( b 'l 

[D--COSP, J 

+ f D-— cosa s ] -21 D-— cos(3 m If D-— cosa, |cos0„ 

v 2 S J l 2 2 S J 

3/2 


, un 7t dvd(3d0 +dvdfd6 

dvd P m = mX o 


_ 7t yl dvd(3 m + dvdff 
Vs ~ 2 


b 

xw, = — cos a, 

s 2 
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The brute-force aspects of the above scheme are obvious. The wake age from 9 = 0 to 180 
degrees is divided into 18,000 segments or 0.01 of a degree. The spanwise segmentation is a little 
more rational. Ninety (90) vortices are trailed, and induced velocity at the wing is calculated in 
between each pair of trailed vortices. This density is far from practical for the real problems 
rotorcraft engineers are solving everyday — with today’s computers. Calculating area as the sum of 
rectangular slivers is hardly advanced. Despite the obvious improvements that might be made, the 
scheme worked. 


The first numerical integration result, compared to the virtually exact solution obtained by 
following Glauert, is shown with Fig. 8. The agreement over 99.9% of the wingspan is more than 
acceptable for engineering purposes. The only problem that occurred was calculation of induced 
velocity at the most outboard span station (i.e., s + 1/2 = 89.5 or xw = 0.4999238476). This span 
station is half way between the vortex trailed from the tip (i.e., m = 90) and the next vortex inboard 
at m = 89. The numerical solution gave v = 19.65 ft/sec versus the virtually exact 1.206 ft/sec. In 
contrast, the port wingtip encountered no such problems. Other than this one ridiculous answer at the 
starboard wingtip (plus being rather slow), the numerical integration scheme functioned in a 
satisfactory manner. 


There are several interesting features to this U-turning wing problem. For example, the 
numerical integration proceeded from port to starboard wingtip with the numerous trailed vortices in 
between. Figure 9 illustrates what the induced-velocity wake age summation appears like for the 
calculation point next to the wing mid-span. The span station of interest is s + Vi (a s = 90.5 degrees, 
xw = 0.0087). Figure 9 shows that dv/dp sees the impending discontinuity at the span station point, 
but the summation averages the calculation to the left of the point with the calculation to the right of 
the point. Therefore, the induced velocity, the integral found from the 

dv 

induced velocity, v = area under — versus 3 

dp 


is obtained with quite reasonable results. 


Another interesting behavior of this U-tuming wing problem is shown with Fig. 10. This 
figure looks at the wake age integration involved with 


— dP = area under 

dp 


dv 

dp 


dp 


dO 


versus 0 


The summation is illustrated at the two vortex trailed wing stations of p = 160 and P = 162 degrees 
(i.e., m = 80 and m = 81). The station at which induced velocity is sought is a s = 161 degrees. The 
fascinating point made here is that the two vortices appear to make enonnous contributions to 
induced velocity within the first 2 degrees of wake age. This is one of the reasons such small 
increments in wake age (i.e., 180 degrees of wake divided into 18,000 segments) are required for the 
rudimentary rectangular area summation scheme. The spanwise distribution of dv/dp, provided with 
Fig. 11, shows just how large the numbers are that, when summed, however, come out v = 
1.1968318858. 
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The situation near the port wingtip, which is closest to the pylon, is quite similar to that near 
the starboard wingtip. The conclusion is that the first 15 degrees of wake age must be very densely 
populated with points if a rudimentary numerical integration is used. Fortunately, advanced methods 
currently used in the rotorcraft industry accomplish the numerical integration much more efficiently. 

Rotorcraft engineers have improved their non-unifonn induced-velocity calculating methods 
by dividing the wake into a near wake and a far wake. Thus Fig. 8 shows that the near wake (i.e., 9 = 
0 to -15 degrees in this case) gives the most trouble for numerical integration, particularly at the 
wingtips. The far wake (i.e., 9 from -15 to 9 = -189 degrees in this case) responds to numerical 
integration with virtually no problems. In fact, the wake age can be divided into 1 -degree segments 
versus 9.91 -degree segments, and virtually the same contribution of the far wake to the induced 
velocity at the wing will be obtained. 


The U-turning wing problem can be divided into a near wake and a far wake with relative 
ease. The induced velocity at the wing due to 189 degrees of circular arc wake in elliptic integral 
form is 


/»+b/2 


(9.43) 


V = 

xw 


dr 

dxv 


4n 


J -b/2 


/’Tt/2 


+ Jo Vl-k 2 sin 2 (^ 


dc) 


XV -xw 


2D + xv + xw 


dxv 


Keep in mind that cj) = tx/ 2 is closest to the wing. Conversely, (j) = 9 corresponds to the end of the 
wake, or where the wing was before it started the U-tum. Therefore, the wake integration amounts to 
integrating from zero up to the start of the near wake, say ())nw, and then adding the near-wake 
contribution, which extends from cj) = cj) N w to cj) = rc/2. It is, however, much more direct in this case, 

to subtract from the total wake (i.e., J q f (4>)dcj) ) the wake behind the near wake (i.e., f (cj))dct) ). 

Thus, the two elliptic integrals are rearranged as follows: 


Near wake = 


(9.44) 


J o -y/l — k 2 sin 2 <|) dcj) - J ( * N " ^/l — k 2 sin 2 <|) dcj) = J * ^/l-k 2 sin 2 <|) d<j> - k 2 sin <|) NM sin i 


Near wake = 


^/l-k 2 sin 2 <|) 


dcj) — 


’4nw 


dcj) = 


\l 1-k 2 sin 2 <j) Jo -\/l-k 


2 • 2 

- sin < 


^ dcj) 


where the angle, x , a new upper limit of integration, is given as 


(9.45) x = 2arctan 


>/ 1 — k 2 sin 2 cj) NW - sin cj) NW Vl-k 2 


coscj^ 
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Then, in elliptic integral shorthand, the near-wake contribution to induced velocity at the wing is 
simply 


( 0 - 46 ) ( V xw)Near = 


Wake 


■+b/2 dr 

dxv 


4n 


E(x> k) - k 2 sin (|) NM sin % | F(x,k) 


XV -xw 


2D + xv + xw 


dxv 


-b/2 


where 


E(x, k) = £ ^1-k 2 sin 2 4> d<() 

1 


(0.47) rx 

F(x,k)= —r , . , 
Jo y]l-\csm-§ 


d(j)v 


In like manner, the induced velocity at the wing due to the far wake becomes 


(0.48) (v„) F „ = 


Wake 


/»+b/2 

- 

'"I’nw 1 

7 = = # 

o yf 1-k sin" cj) 

8T 

dxv 

£ nw -Jl-k 2 sin 2 cj) d(j) 

H - 

4n 

XV -xw 

2D + xv + xw 

J -b/2 




dxv 


where, to repeat, k 2 = 


4(D + xw)(D + xv) 
(2D + xv + xw)“ 


This completes the discussion of the U-turning wing problem. Two points have been made. 

(1) A virtually exact calculation of induced velocity without numerical integration questions has 
been provided. Figure 6 gives an example to which any numerical integration result may be 
compared. 

(2) The so-called near wake is built upon the fundamental integral solved by the fixed-wing 
community, which is 


(0.49) 


Fixed Wing (v xw ) 


C+bl2 ST 


oxv 

1 

4n 

J -b/2 

xv -xw 


dxv 
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II. The Wing Done With Rotor Notation 

The next step along the bridge from the fixed-wing world to the rotary-wing world is to 
change to rotor blade geometry and notations. This geometry is provided with Fig. 12. The intent 
here is to think of a one-bladed rotor in hovering flight. (The case of forward flight begins in Part III 
of this note.) The rotating wing is again represented as a lifting line having an elliptical bound 
circulation. The trailing vortex wake structure resembles a lock washer or, perhaps more 
descriptively, a “slinky” spring-like toy. That is, each vortex will have a circular path just like the U- 
turning wing problem. The addition is that each trailing vortex spirals downward at a constant rate 
and the wake age can extend back to the beginning of time. In short, the wake of this rotating wing 
or, better yet, a one-bladed hovering rotor, need not stop at minus 180 degrees as was done with the 
U-turning wing. 

Now study Figure 12’s geometry quite closely. The rotor blade is shown rotating around the 
Z axis, which is positive down. The blade rotates in the X - Y plane where X is taken positive 
“forward,” although what forward means in this case of no forward speed is not really meaningful. 
The Y axis is normal to the X - Z plane. This X, Y, Z axis system does not rotate with the blade, nor 
does it move with time. 

The blade itself has a radius, R, measured from the Z axis. The blade is shown in Fig. 12 
with a root cutout, r c . Thus, the blade’s span (thinking in terms of a wing) is simply R - r c . Any 
radial station, measured positive outward from the Z axis and along the lifting line, is denoted by, r. 

The blade’s rotation angle is measured by the azimuth angle \|/, which, for convenience, 
equals zero when the blade lies along the negative X axis. The azimuth angle increases with time 
simply as the rotational speed Q in radians per second times time in seconds. Obviously, Fig. 12 is a 
snapshot at any given time you care to start the watch. The blade is simply going round and round 
and the trailed vortex structure is left in space to descend at a prescribed rate, dZ/dt, in feet per 
second. Figure 12 is drawn to imply that the blade has been rotating for the time it takes to complete 
about two revolutions. However, the arrow head shown at the end of the one trailed vortex implies 
that time has been going on forever and, with a longer piece of paper, the spiral would extend down 
and around the Z axis to infinity. 

As with the preceding fixed-wing examples, the Biot-Savart law requires a very careful 
mathematical definition for the blade and vortex wake dimensions. This is even truer when tackling 
the rotor blade problem. Using Fig. 12 as the geometric model of a single-bladed rotor, the radius 
station, r, at which the induced velocity is sought is set relative to the X - Y - Z axis by 

(0.50) Xr = -rcos\|/ Yr = +rsin\|/ Zr = 0 

A vortex trails aft from the blade from radius station rv. Any segment of this long spiraling vortex is 
therefore deposited at a point behind the blade and below the X - Y plane. This vortex segment, dS, 
is located at the coordinates 
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1 ry 

(0.51) Xv = -rvcos0 Yv = +rvsin9 Zv = H -At 

dt 

Equation (0.51) introduces the possibility that the wake trailing behind the rotor blade does not stay 
in the X - Y plane in which the rotor blade turns. This is quite different from classical fixed-wing 
wake geometry assumptions. The rotor-blade wake descends with some velocity, dZv/dt, and this 
velocity need not be constant. In fact, in the more complete analyses of rotor systems, freedom is 
given for any given vortex segment to wander throughout the X - Y - Z axis system. For elementary 
discussion purposes, just assume a constant-diameter spiral (i.e., neither rv nor dZv/dt are influenced 
by time). Furthermore, it is not necessary at this point in the discussion to be more specific about the 
vortex segment’s descent velocity. 


Now as to the matter of time introduced by Eq. (0.51); for the rotor blade, the time 
increment, At, represents the time, t, it takes to travel back from the blade, which is at the snapshot 
azimuth angle, \\i, (at time, t 0 ) to the vortex segment, dS, which is located at 9 (at time, t) . That is 

(0.52) At = t - 1 0 


Since the blade’s angular rotation speed is Q, it follows that 
9 -\| / 


(0.53) At = t - 1 0 = - 


Q 


which immediately says that 

r -7 dzv . dzvf9-\i/^ 

(0.54) Zv = + At = - 

dt dt l Q 


dzv 9 _ dzv Q 


Qdt 


dO 


In many propeller studies, dZv/d0 is a measure of helix angle or sometimes propeller pitch. 


Now, quite methodically, begin with the Biot-Savart law as previously stated: 


(0.55) d(dv x /dxv) = — 

47t 


LjdSj -LjdS; 
L 3 


d0 


where the dimensions are expressed as vectors. The distance, L, between the vortex segment to the 
radius station, r, at which the induced velocity is sought, is determined by 

L = Lj + Lj + L k = (Xr - Xv)i + (Yr - Yv) j + (Zr - Zv)k 
where Xr = -rcos\j/ Yr = +rsin\j/ Zr = 0 

Xv = -rvcos9 Yv = +rvsin9 Zv = -T^^0 

d0 

L 3 = [(Xr - Xv) 2 + (Yr - Yv) 3 + (Zr - Zv) 2 J' 2 


(0.56) 


and 


In a similar manner, the vortex is described by 


18 



(0.57) S = S,+S J +S k =(Xv)i + (Yv)j + (Zv)k 

but with rotor-blade vortex geometry, 

dS; = ^-d0 = ^^d0 = (rvsin0)d0 
1 50 50 v ’ 

5S ■ 5 Y v / \ 

(0.58) dSj = — 1 d0 = - — d0 = (rvcos0)d0 

J 50 50 v ’ 

ds„ =-^L d e=— e 

50 50 

The substitution of this rotor-blade geometry into the Biot-Savart law gives, with simplification 


(0.59) d(dv r /drv) = — 

4n 


rv 2 -r(rv)cos(\|/-0) 


rv 2 + r 2 - 2r (rv) cos (\|/ - 0) + 


5 dzv^ 


Vd9y 


3/2 


d9 


Because this is a single -bladed rotor in hover, there is an opportunity to choose the blade 
azimuth position at any value that is convenient. This is because the hover problem is completely 
symmetrical around the Z axis and “forward” has no meaning. (This is not possible in the forward 
flight case, as you will see in Part III of this note.) Given this latitude, choose \\i = 0. This reduces 
Eq. (0.59) to 


(0.60) d(dv r /drv) = — 

47T 


rv' 


-r(rv)cos0 


rv 2 + r 2 - 2r (rv) cos 0 + 


f dzv^ 


V dOy 


-| 3/2 


d0 


which is a variation on Eq. (0.29) arrived at for the fixed wing after a U-turn. The variation is, of 
course, the addition of the term 


2 dzv^| 


0 2 


Now let me address the wake age parameter, 0. The vortex wake leaves the blade lifting line 
referenced to \j/ = 0. Therefore the smallest that 0 can be is 0 = \\> = 0. The vortex spirals backwards 
(round and round and down) to 0 = -oo. While it is not practical to go all the way back to 9 = -oo, the 
intent is to go as far back as possible to avoid missing any influence of what many refer to as the far 
wake. Nevertheless, the integration over the wake age becomes 
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(0.61) dv r /drv = 


4n 


rv 2 -r(rv)cos0 


rv 2 + r " - 2r (rv) cos 0 + 


^ dzv v 


V d0 y 


— 13/2 


•dO 


The immediate objective is to perform the integration required by Eq. (0.61). Suppose the 
vortex circulation, y v , does not vary with wake age. And suppose the vortex descent measure, 
dzv/d0, is constant. The only immediate numerical problem that is apparent in Eq. (0.61) is when the 
integrand’s denominator is identically zero or so close to zero that numerical integration built into 
Mathcad flounders. This situation will only occur when 0 is actually zero and when r = rv. 

To illustrate the results of integrating Eq. (0.61) over the wake age using Mathcad’s built-in 
scheme, choose some rational values for a single-bladed rotor. For this example 


Parameter 

Value 

Rotor radius, R, feet 

30.0 

Root cutout, r c , feet 

3.0 

Trail vortex from radius, rv, feet 

27.0 

Calculate induced velocity at radius, r, feet 

16.5 

Trailed vortex strength, y v , ft 2 per second 

1,265.0 

Vortex descent measure, dzv/d0, feet per radian 

2.0 

Wake age, 0 en d, radians 

-157.0 

Number of 360-degree spirals 

25 

Distance from rotor down to last spiral, zv, feet 

314 


Mathcad was quite happy to whip out the answer so long as the input was chosen that avoids r = rv 
by a considerable margin. (Just imagine the wake age integral in summation form with wake age 
segments of 0.01 degrees! It would take forever.) In the blink of an eye Mathcad produced Fig. 13. 
Figure 13 shows that induced velocity accumulates very quickly with increasing number of spirals 
included in the integration. Closer inspection of Fig. 14, an enlargement of Fig. 13, shows that the 
vortex segments included in the first half spiral contribute in a rather linear fashion to the induced 
velocity at blade radius station, r. But as the spiral continues around, it comes back underneath the 
blade (see Fig. 12) and the close proximity allows the vortex to become very influential. In short, the 
denominator in Eq. (0.61) with 0 around -270 degrees is very influenced by the vortex descent 
measure, dzv/d0. 

In advanced, modern numerical integration of the real rotor system problem, wandering 
vortices do come close and, indeed, even intersect the lifting line or lifting surface of the blade that 
trailed the vortices. That is to say Eq. (0.61)’s denominator finds some way to go to zero in the real 
problem. Furthermore, when there are more blades to consider, the ensuing possibilities of this 
denominator going to zero has been a major source of frustration in achieving accurate solutions for 
the rotor hovering and forward flight problems. The current crutch most widely used has been to add 
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what is called a vortex core to the denominator of Eq. (0.61) so zero can never occur. The fact that a 
vortex does have a real physical core of measurable diameter makes the core’s inclusion in Eq. 
(0.61) reasonable. This concept of a vortex core diameter, VCD, means that Eq. (0.61) can be 
written as 


(0.62) dv r / drv = 


4n 


rv 2 -r(rv)cos0 


rv 2 +r 2 -2r(rv)cos9 + 


f dzv^ 2 


-i3/2 


*d0 


0 2 +VCD 2 


To make progress over recent years, technologists have achieved considerable numerical stability 
even when tackling the hardest problems using this concept of a vortex core. For this part of the 
note, completely dismiss the whole vortex core diameter issue. 


Going on then, consider the collection of all the trailed vortices by performing the radial 
integration over the blade span from the root cutout, r c , to the blade tip, r = R. This step begins by 
assuming a bound circulation distribution carried on the lifting line. Assume the bound circulation 
has an elliptical distribution 4 described by 


(0.63) 



(2r-<;-Rf 

(r-A 


which gives the circulation strength of a vortex trailing from radius station, rv, as 


(0.64) 


Yv 


-2RT o 

2rv - r c - R 

, (2rv-r c -R) 2 

. ( r -A 

1 ( R -A 



>drv 


The induced velocity along the blade’s radius collecting all vortices is then 


(0.65) 


2 r 


4ti 


1- 


2rv-r -R 


(2rv-r c -R) 2 [_ (R~r c ) J 

(R-O 2 


rv 2 -r(rv)cos0 


id0drv 


rv 2 +r 2 -2r(rv)cos0+ 9 


Now, let me illustrate a completely numerical solution for this double integral. 


4 Keep in mind that the more general distribution used by fixed-wing engineers, Eq. (0.23), could be used to extend the 
results beyond an elliptical distribution. 
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Now, from Parts I and II of this note, the numerical integration of Eq. (0.65) will not 
converge to a solution because of any number of singular points. Mathcad’s built-in integration 
scheme flounders even if a close-proximity situation occurs. You also know that the induced 
velocity — given an elliptical bound circulation along the lifting line — must contain the solution for a 
fixed wing. That is, at the very least 


v, = ■ 


The approach then is to first borrow the fixed-wing solution technique of letting 


r = 


R + r 


R-r, 


(0.66) rv = 


drv = 


f R+O 


fR-O 

l 2 J 


l 2 J 


R-r 


cos a = a-bcosa 
cosp = a-bcosp 


sin p dp = bsinpdp 


where P goes from 0 to n. This coordinate system change at least transforms the vortex circulation 
strength, y v , of Eq. (0.64) to the very much simpler 

(0.67) y v =-r o cospdp 

The completed substitution of Eqs. (0.66) and (0.67) into Eq. (0.65) restates the induced-velocity 
double integral problem as 

( 0 . 68 ) 


r n r° 

— r 0 cosp 

(a -bcosa)(a-bcosP)cos9-(a-bcosP)~ 



4 71 

J 0 j -CO 

(a-bcosP)~ + (a -b cosap -2(a-bcosa)(a-bcosP)cos9 + 

f dzv^ 

IdeO 

2 1 

02 

3/2 


So, now to numerically integrate Eq. (0.68). The scheme allows Mathcad to perfonn the 
wake age integral, but the spanwise integration is perfonned with a summation. This takes advantage 
of Mathcad’s very, very fast integration over wake age. Thus, a useable numerical integration 
scheme reads like this: 
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INPUT 


r =2 R = 1.2 r = 0.2 — = 0 

0 c d0 

Dimension integration. 

Extent of wake age, WA = -n 

Number of radial stations at which vortices leave rotor blade, M = 90 
Range of radial stations where induced velocity is calculated, s = 0, 1 M - 1 


Then proceed with these calculations 


R-r 


b = 


R + r 


n ( 1 

a = — s + — 
Ml 2 


m = 0,1 M 


Pm = m 

Pm 1 M ' 


dvdp m = 


-r o cosp 


4rt 


(a -b cos a) (a -bcosp)cos9-(a -bcospf 


(a-bcosp)' + (a-bcosa)~ -2(a -bcosa)(a -bcosPlcosO + j^ 1 ^ j 0 2 


>d© 


y 7t y dvdp„, +dvdp„ 
s Mti 2 
r = a - b cos a. 


Notice on this sample input that it corresponds to the U-tuming wing problem of D = 0.7, b = 1, r o = 
2, which gives the induced velocity distribution shown on Fig. 8. This check case did check. 
Mathcad was faster and more accurate with the wake age integral broken into several ranges. The 
following seemed best 

C f («-P)de=f JLi f(a,p)de + J^ i f(a,p)de + Jj(a,p)de+J w y(a,p)de 

180 


Now consider a practical case of a 4-bladed, 44-foot-diameter rotor lifting about 10,800 
pounds while operating at a tip speed, V t = OR = 600 ft/sec. The above input requires an estimate for 
maximum bound circulation, T 0 , and the vortex segment descent rate, dzv/d9. Rational estimates for 
both parameters can be obtained. To begin with, assume each blade has a root cutout of r c = R/6 and 
each blade has an elliptical bound circulation over its span. The bound circulation’s maximum will 
be T 0 = 225 ft 2 /sec, which is calculated from 


7t 


(0.69) Lift per blade = — pRV t T o 

8 


1- 


f rc^~ 


vRy 


so 


r = 


8L 


7tpRV t l-(rc/R) 
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Next, assume a vortex segment, dS, is carried downward with the downwash velocity. From simple 
momentum theory, this downwash or induced velocity is 19.3 ft/sec when calculated with 


(0.70) Momentum induced velocity = v ; = 


L 


' 2p7tR" 


Then, in one second, a segment will have descended 19.3 feet. But in one second, the rotor blade 
will move through a wake age angle of 9 = Q times one second. Therefore, the vortex segment 
descent rate is simply 


(0.71) 


dzv V; times 1 second 


V; 


d9 Q times 1 second 


— ^ — = R — 

V,/R v, 


which is roughly 9.7 feet per radian for this representative example. This means that in one complete 
spiral revolution, the vertical distance will be 0.7(2 tt) feet or roughly 4.4 feet for every revolution. 
This is a measure of Fig. 12’s spiral spacing or pitch. 


To summarize this typical one-bladed rotor’s operating situation, the representative 
calculations that follow were based on 


INPUT 

r =225ft 2 /sec R = 22ft r. =R/6 — = 0.7 ft /rad 

0 c de 

Dimension integration. 

Extent of wake age, WA = variable for this example 

Number of radial stations at which vortices leave rotor blade, M = 90 

Range of radial stations where induced velocity is calculated, s = 0, 1 M - 1 


An immediate question is: How many spirals does it take to accurately approximate an 
infinite wake? As seen from Fig. 15, perhaps about one hundred, 369-degree spirals captures the 
problem in a promising way. This amounts to 449 feet of wake or about 19 rotor diameters. A more 
quantitative measure of a practical engineering solution is induced horsepower, HPj. This major 
contributor to total power is calculated as 

(0.72) HP, =-t-f R v r dL r = -I-| R v I (pV r r,)dr 

The calculation was made using the summation of 99 rectangular slivers. The results, tabulated 
below, indicate that even 29 spirals will neglect about 2 to 3 horsepower, which is on the order of 2 
to 3 percent missing horsepower because the wake age was not extended to infinity. The ratio of 
calculated induced horsepower to ideal momentum horsepower (i.e., Lvi/559 = 95.524 horsepower) 
is the fourth column in the table. The reciprocal of this ratio, known as Figure of Merit, is provided 
by the sixth column. A Figure of Merit above 1.9 is not possible and thus at least 19.5 spirals are 
required before a rational answer starts to become apparent. Clearly, an elliptical bound circulation, 
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which is ideal for a fixed wing (i.e., the equivalent of FM = 1 for a rotating wing), is far from ideal 
for a rotating wing. 


Number 




Ideal 



of 

Blade 

Calculated 

Induced 

Calculated 

Figure of 

Spirals 

Lift 

Induced HP 

HP 

HP/HPj deal 

Merit 

0.5 

2712.52 

29.57 

95.52 

0.31 

3.23 

1.5 



55.31 



0.58 

1.73 

2.5 



70.41 



0.74 

1.36 

3.5 



80.32 



0.8408 

1.189 

5.5 



92.33 



0.9666 

1.0346 

10.5 



104.44 



1.0934 

0.9146 

20.5 



110.25 



1.1541 

0.8665 

100.5 



112.79 



1.1808 

0.8469 

200.5 

1 

< 

112.84 


1 

1.1813 

0.8465 


Do not assume that this example has provided anything more than a crude estimate for a single- 
bladed rotor’s actual performance. The prescribed wake geometry used is very, very far from the 
wake visually observed in any number of experiments. The example would surely come out 
differently — and more accurately compared to experiment — if (or when) this sample problem is 
calculated with advanced methods in use today. 

For the sake of completeness, Fig. 16 gives the running lift, AL r , and running horsepower, 
AHPi, versus radius station. These data are for 100.5 spirals. 

This concludes Part I and Part II of this note. The next step is to apply this background to the 
single-bladed rotor in forward flight. 


25 



26 



III. A Rotating Blade in Forward Flight 


The next step in the bridge between a fixed and rotating wing is examined with this portion 
of the note. The rotating wing in forward flight deposits a rather complicated vortex wake in space 
as it flies away. A hint of the wake’s complexity is seen below. 
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In the sketch (pg. 27), the rotating wing — the single-bladed rotor — is rotating about the 
Z axis in the X - Y plane at angular velocity, Q , in radians per second. The blade is attached to a 
hub. The hub is located at X=Y=Z=0. The X - Y - Z axis system moves through space straight 
along the plus X direction with forward velocity, V, in feet per second. The axis system neither 
pitches nor rolls and it does not climb or descend. The rotor blade has a radius, R, in feet and a root 
cutout, r c , in feet. The tip speed of the rotor is V t = OR in feet per second. The one vortex shown 
illustrates the drifting-down, spiraling path typical of all vortices. 

A more complete picture of the wake complexity is provided with Figs. 17 through 20. These 
figures have been drawn for an advance ratio, p = V/V t = 0.5, and show the planform view. Start 
with Fig. 17 where the blade is at the azimuth angle, \| /, of 90 degrees. This azimuth is generally 
referred to as the advancing side of the rotor disc. In the fixed-wing problem, all vortices trail 
straight back, parallel to the X axis. Obviously, this is not true for a rotor blade. The tip vortex traces 
out a prolate cycloid, while the root cutout vortex follows a curtate cycloid. Figure 18 shows the 
planform view with the blade at the 180-degree azimuth position. Figure 19 places the blade at \(/ = 
270 degrees, which is the retreating side of the rotor disc. Notice that the two outboard vortices 
trailed from the blade sharply turn nearly back on the blade itself. Finally, Fig. 20 shows the blade at 
its most downwind position of \j/ = 360 or zero degrees. Here the trailed vortices really attack the 
blade. The possibilities of any given vortex directly intersecting the generating lifting line are quite 
real in the practical problem. 

Before bringing the Biot-Savart law to bear on this problem, there are a few aspects of the 
notations to observe. To begin with, rather than deal with the geometry dimensionally, use the 
conventional rotor non-dimensional notations of 


x = — radius station where induced velocity is sought 


(0.73) 


xv = — radius station where vortex is trailed 
R 

a t = angle of attack of the tip path plane, positive nose up 


V 

p = — cos a 

V 


tpp 


advance ratio 


V; = induced velocity calculated by momentum theory 

Vsina ( „„ — v ; 


^tp P 


tpp 


V 


— inflow ratio 


Now, using Fig. 21 as the reference, methodically begin with the Biot-Savart law as 
previously stated: 


(0.74) d(dv x /dxv) = — 

4n 


LjdSj -LjdS, 
L 3 


d9 
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where the dimensions are expressed as vectors as displayed on Fig. 21. For the distance, L, between 
the vortex segment, dS, to the radius station, x, at which the induced velocity is sought, you have 


L = L; + L . + L k = (Xr - Xv)i + (Yr - Yv) j + (Zr - Zv)k 
where Xr = -R[xcosv|/] Yr = +R[xsinv|/] 

Xv = -R[xvcos0 + p(v|/-0)] Yv = +R[xvsin0] 

and L 3 = [(Xr - Xv) 2 + (Yr - Yv) 3 + (Zr - Zv) 2 ] 3/2 


Zr = 0 

Zv = -RZ tpp (v|/-0) 


In a similar manner, the vortex is described by 
(0.76) S = S,+S j +S k =(Xv)i + (Yv)j + (Zv)k 
but with rotor-blade vortex geometry, you have 

PC Pv v 

dSj = — d0 = — — d0 = +R(xvsin0 + p0)d0 
1 50 50 v ’ 

5S • 5Yv / x 

(0.77) dS ; =^d0 = d0 = +R(xvcos0)d0 

J 50 50 v ’ 

5S 

dS k = — ^ d0 = +R \ 0 

k 50 pp 


The substitution of this rotor-blade and wake geometry into the Biot-Savart law gives, with 
simplification 


(0.78) 

LjdSj -LjdS; 

V 


xv 2 -x(xv)cos(\|/-0) + p[xvsin0-xsin\|/ + xv(\|/-0)cos0] 
r|[xv 2 + x 2 -2x(xv)cos(\|/-0)] + 2p(xvcos0-xcos\|/)(\|/-0) + (\|/-0) 2 (p 2 +Z( pp ) 


As an intermediate step, integrate the Biot-Savart with respect to wake age assuming the 
vortex circulation is simply a constant. That is, tackle 


(0.79) (d\ 


Yv. 

c V 

1 

r 

"a. 

JZ) 

1 

r 

'~£u 

yi 

i 

47T 

«/ 

-00 

L 3 


dO 


The purpose of this step was to test Mathcad’s built-in integrator and be sure it did not flounder at 
any azimuth, \| / , or radius station, x, with the provision that x ^ xv. To perform this test, I set y v = 
4k, and R = 1. I found that the near wake needed integration in several parts and finally selected 
integration as follows: 
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rv 


(0.80) 


LjdSj - LjdS; 

L? 


7^1 71 _ ^ < c 

r vi/ /•viz 1 /•vi; 2 /•viz 5 (• vi/ 15 

d9 = J +j ” +| ” + j ”, +| ” + 

J vj/ 1 J vi/ 2 J vi/ 5 J vi/ 15 J vi/ 90 

180 180 

py-Ti p- 

L + f 


1 


ISO 


v|/ 90 


•I vi/ 15 

180 


)w ~— 90 
180 


u / 180 

180 


-571 

.1-2077 


Mathcad’s built-in integrator struggled with the calculation at and near the blade’s root end (i.e., 
around the root cutout region) in the azimuth region from vj / = 330 to 360 degrees. The reason for 
this struggle is, of course, the near-zero value of L, as Fig. 20 clearly shows. Adding a vortex core 
diameter would obviate the problem in regions where the vortices are so closely packed. 

The next step requires picking a bound circulation for the blade’s lifting line. Suppose the 
elliptical distribution is chosen and this distribution does not vary with azimuth. This means there 
will only be trailed vortices and no shed vortices to add to the problem for this example. (The case of 
a shed wake will be addressed shortly.) Thus, 


(0.81) 


r =r 

X,\\l 


1 


, J'2x-x, -I f 

(1-A 2 


does not vary with azimuth, \j / 


Then, borrowing from the fixed-wing solution technique, let 

I cos a = a-bcosa 


x = — = 

R 


fi+o 


f i-O 

l 2 J 


l 2 J 


(0.82) 


xv = — = 

R 


fl + x ) 



l 2 J 


l 2 J 


cosP = a-bcosp 


dxv = 


1 — x„ 


sin p dp = bsinpdp 


where p goes from 0 to n. This coordinate-system change at least transfonns the trailed vortex 
circulation strength, y v , to the very much simpler 

(0.83) Y v = -T 0 cos P dp 

and thus the induced-velocity double integral problem for an elliptical bound circulation that does 
not vary with azimuth is restated as 


(0.84) 


v„ = 


r n 



-r 0 cosp 


4tcR 

J 0 

-00 


xv 2 - x (xv)cos (\|/ - 9) + p [xvsin 0 - x sin \\i + xv(\|/ - 9) cos 0] 


{[xv 2 + x 2 -2x(xv)cos(\|/-0) +2|i(xvcos9-xcos\|/)(\|/-0) + (\|/-0p (p 2 +A,jj, p ) 


d0dp 
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The question now arises as to what value of the maximum bound circulation, r o , is 
representative for this example. The calculation of this single-bladed rotor’s average or steady lift 
follows as: 


dL rv = p(V rH , xf r v )dr both V and T are vectors 

dL x , v = PV X / S , ¥ Rdx = p[V t (x + psin\(/)] 


r 


1 


(2x-x c -l)~ 

(l-x c ) 2 


Rdx 


r 1 


L ,= 


(0.85) 


j >*,= 


p[V,(x + psin\)/)] 


r 


1 


1 ( 2 x ~ x c -!) 2 

( 1 — c ) 2 


Rdx 


L v =|p R V t r 0 (l-x c )(l + 2qsinv + x c ) 


1 r 2 n 1 C 2n n , 

L steady=^f 0 L v d ¥ = — - P R V t T 0 (1 - X c ) ( 1 + 2p sin \|/ + X £ ) d V 

Zk Ztz Jo o 


8L 


Steady 


Steady Lift per blade = — pRVX [l-x^l so T o = r 

8 P ° L cJ 7rpRV t [l-x;] 


On this basis, choose the blade geometry from Part II 


R = 22 feet 

x c = 1/6 non-dimensional 


and the forward flight conditions for this Part III problem are: 


p = 0.002378 slugs per cubic feet 
V t = 603.605 feet per second 
p = 0.5 non-dimensional 
^-tpp = - 0.03 non-dimensional 
T 0 = 225 square feet per second 

in which case the steady lift per blade is 2,712 pounds at 178 knots. 

With this information as input, Mathcad and its built-in integrator is used to calculate Eq. 
(0.84) as follows: 
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INPUT 

11 = 0.5 T o = 225 R = 22 x c =l/6 L tpp =-0.03 

Dimension integration. 

Select azimuth station, v[/ = variable for this problem 
Extent of wake age, WA = -20rt 

Number of radial stations at which vortices leave rotor blade, M = 90 
Range of radial stations where induced velocity is calculated, s = 0, 1 M - 1 


Then proceed with these calculations 


1-x. 


a = — s + — 
M 


Xs R 


1 _ 

2 

1 + x. 


b = 


l-x„ 


l + x. 


cosa, = a-bcosa„ 


m = 0,1 M 

l+x, 

xv m = — = ' 

m R 


P,„ = — m 
Pm 1 M ' 


1-x, 


dvd[T m = 


-r 0 cosp 

4tiR 


c°sp m =a-bcosp m 

xv 2 - x (xv) cos (\|/ - 9) + p [xv sin 9 - x sin \|/ + xv(\|/ - 9) cos 9] 


| xv 2 +x 2 - 2x(xv)cos(\|/-9)] + 2p(xvcos9-xcos\|/)(\|/-9) + (\|/-9)‘ (p 2 + Lj pp )j 2 


d9 


n yl dvdp s m + dvdP s 
^ 2 
x = = a -bcosop 


The calculation of induced velocity at the lifting line — as created by the trailed vortices from 
the lifting line — now follows from the above scheme. Additionally, lift and horsepower distributions 
are calculated from 


( 0 . 86 ) 


AL s , ¥ =[pV t (x s +psint|/)R]<r jl- ( 2 * s Xc ^ 

V ( 1 _X c) 

M-l 

L ¥ =X AL s => Single-blade lift at input azimuth 

s=0 

AHP S . V = ^ v s V AL S ¥ => Running Horsepower 

M-l 

HP ¥ = ^ AHP S ¥ => Single-blade Horsepower at input azimuth 

s=0 

where A , = I 


[ A s => Running Lift 


fl + O 


f n ) 

( n , .3 


cos 


-cos — s + 1 

l 2 ) 


Im J 
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A representative illustration of azimuth varying induced velocity at several radial stations is 
shown on Figs. 22 and 23. Remember that this result is for a single-bladed rotor lifting 2,712 pounds 
at 0.5 advance ratio (i.e., V = 178 knots at V t = 603.6 fps). Even more interesting is the azimuthal 
variation of the blade’s total lift and induced horsepower, shown on Fig. 24. The average or steady 
induced horsepower is obtained by 

(0.87) Steady Induced Horsepower = HP; HP 1|( 


The result is HP; = 50.1 hp. This induced horsepower, calculated with the prescribed wake, 
compares to the ideal induced horsepower calculated in Part I as 


( 0 . 88 ) 


Ideal Induced Horsepower = Ideal HP ; 



' Lift ^ 
v 2p7tR 2 V y 


for jlx > 0.2 


which yields ideal HP; = 6.13 hp. This means that the elliptically loaded, single-bladed rotor requires 
about 8 times the power calculated by simple momentum theory! 


An additional point made by Fig. 24 is that the highest lift is carried primarily on the 
advancing side of the rotor disc. As such, the rotor is out of trim because of the rolling moment 
inferred by Fig. 24. 


This leads to the question: If the rotor has an elliptical bound circulation that varies with 
azimuth so that the rolling moment is zero, what is the induced power? And this leads to the fourth 
and final Part of this note. 
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IV. A Rotating Blade in Forward Flight With Zero Rolling Moment 


The last step in the bridge between a fixed and rotating wing is concluded with this portion 
of the note. A fixed wing generally flies with zero rolling moment, RM. This same criterion can be 
applied to a rotating wing. [I refer to a rotating wing with RM = 0 as a “balanced” rotor.] Suppose 
the bound circulation of the single-bladed rotor is described by 


(0.89) r i¥ =(r 0 + r,sinv|i) 

In this case the lift and rolling moment are found as 


i (2x-x c -l) 
(1-xJ 2 


(0.90) 


L st ea dy = £ pRV t r o [1 - X c ] [(i + x c ) r o + pr, ] 


RM 


Steady 


71 

128 


pR 2 V t r o [l-x c ][8|a(l + x c )r o +(5 + 6x c + 5x’)r i " 


and if the rolling moment is set to zero 

8p(l + x c )r o 

(5 + 6x c +5x[) 


r> = - 


(0.91) 


L St ea dy =-pRV t r o [l-X^] 


1- 


8p 2 


(5 + 6x c +5x 2 ) 


r 


The Biot-Savart law of Eq. (0.84) now must include the trailed vortex circulation, y v , which 
is azimuth varying according to Eq. (0.89). Therefore, for this Part IV of the note: 

(0.92) 



r n 

-(r o +Tj sin\|/)cos(3 

xv 2 -x(xv)cos(\|/-0)+ |.i[xvsin0-xsin\|/ + xv(\j/-0)cosO] 

^X,\| 1 

J 

o J 

4jtR 

—00 

3 

{[xv 2 +x 2 -2x(xv)cos(\|/ — 0)] + 2|a(xvcos0-xcos\|/)(\|/ — 0) + (\|/ — 0)~ (|F +A.[, p )| 2 


The numerical integration of this slightly different equation (i.e., with azimuth varying circulation) 
follows exactly that given in Part III of this note. 

The contribution of trailed vorticity to induced velocity at the lifting line of this “balanced” 
rotor is shown in Figs. 25 and 26. The total blade’s lift and induced power is given with Fig. 27. 
These results are for the balanced rotor producing 2,712 pounds of lift at 178 knots and can be 
compared to Fig. 24, which is for the “unbalanced” rotor. Only trailed vortices are contributing to 
the induced power at this point. The additional induced velocity and horsepower due to the shed 
wake are still to be included. 
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A numerical solution scheme for the shed wake is relatively simple compared to that for the 
trailed wake. However, in the shed-wake problem, the vortex left behind the blade has a radial 
geometry and circulation just like the blade’s when the blade was at that azimuth. However, the 
vortex circulation is the negative of the blade’s bound circulation when it was there. The geometry 
of this shed-wake problem is shown on Fig. 28. Again, the conventional rotor non-dimensional 
notations are: 


r 

x = — radius station where induced velocity is sought 
r 

xv = — radius station where vortex is trailed 
R 

a tpp = angle of attack of the tip path plane, positive nose up 

(0.93) y 

u = — cos ol advance ratio 

r y W 

V t 

v ; = induced velocity calculated by momentum theory 
VsinoL -v. 

Z trai = — inflow ratio 

pp V t 


Now, using Fig. 28 as the reference, begin with the Biot-Savart law 


(0.94) d (dv x / dxv) = — 

4n 


LjdSj -LjdSj 

L 3 


d9 


where the dimensions are expressed as vectors as displayed on Fig. 28. The distance, L, between the 
shed-vortex segment, dS, to the radius station, x, at which the induced velocity is sought, is 

(0.95) 

L = L; + L . + L k = (Xr - Xv)i + (Yr - Yv) j + (Zr - Zv)k 

where Xr = -R[xcosi|/] Yr = +R[xsinv|/] Zr = 0 

Xv = -R[xvcos0 + p(v|/-0)] Yv = +R[xvsin 9] Zv = -RA lpp ( ¥ -e) 

and L 3 = [(Xr - Xv) 3 + (Yr - Yv) 3 + (Zr - Zv) 3 ?“ 

3 

so L 3 = R j[xv 2 +x 2 - 2x(xv)cos(v|/ -0)] + 2p(xvcos0-xcosv|/)(v|/-0) + (v|/ -0)" (p 2 +Z( pp )j 2 

Notice that the length, L, is the same for the shed wake as for the trailed wake. To continue then, the 
vortex element is described by 

(0.96) S = S 1 +S J +S k = (Xv)i-r(Yv) j + (Zv)k 

but with rotor-blade shed-vortex geometry, the vortex extends radially and so the partial derivatives 
of S are with respect to xv. Thus 
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(0.97) 



dS 

axv 

dS, 

= — L dxv = 

dxv = -R cos 9 dxv 

l 

dxv 

dxv 


5Sj 

oYv 

dS, 

= — - dxv = 

dxv = +R sin 0 dxv 

J 

dxv 

fixv 


<3S, , 

dZv , 

dS k 

= — -dxv = 

dxv = 0 


dxv 

Sxv 


With this information, the geometric part of the Biot-Savart law becomes 


(0.98) 


i 

r 

~Cu 

_yi 

1 

r 

'~'o- 

yi 

i 

x sin (\j/ - 0) + p (\|/ - 9) sin 9 

L 3 

i j 


dxv 


The Biot-Savart also needs the shed-vortex circulation as it varies with wake age. Thus, 

( dr „ ^ 


(0.99) Shed y v = - 


v % j^r 9 


d9 = -(Tj cosO) 


(2xv-x c -l)- 

0-*,) 2 


de 


and the double integral giving induced velocity at any radius station, x, and azimuth, \( / , using Eqs. 
(0.98) and (0.99) is 


(0.100) v X;V = 

f 1 

r v 

( r iCOS0) 1 

t/ 

x c * 

-00 



(l-x .) 2 


x sin (\|/ - 0 ) + p (\j/ - 0) sin 9 


L 


d0dxv 


It is particularly important to study the integration with respect to wake age, 0, before 
discussing a complete integration of Eq. (0.100). The reason is that this portion of the integration has 
a definite possibility of “blowing up.” Therefore, reverse Eq. (0.100)’s integration order so that it 
reads as 



c V 


* 1 

- 

(0-101) v x>v = 

< 



( r l COS0 )J 

»/ 

-00 

J 

x c 

V 


(2*v-V-lj 
0-*,) 2 


x sin (\j/ - 0) + p(\j/ - 0) sin 0 


L 3 


dxv 


dO 


In the first place, the integration is not a problem when 9 = \j/ because the numerator of the integrand 
is zero (i.e., \|/ - 9 = 0). Furthermore, from Eq. (0.95), the distance, L, between the vortex segment 
and the point on the blade where induced velocity is sought reduces to 


3 

L 3 = R j[xv 2 +x 2 -2x(xv)cos(0)]jr =R j(xv-x)" 
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which is only zero when x = xv. The physical meaning of this situation when 9 = \j/ is that the only 
vortex that exists is the blade’s bound circulation and this straight-line vortex can not induce a 
velocity on itself. 

Now look at the solution when a shed vortex is in the region \|/ - 9 > 0 or, if you prefer, in 
the near wake. Suppose, for example, that the blade is at the 135-degree azimuth and 9 is in the 
range 135 degrees backwards to 132 degrees in 1/39 of a degree increments. This represents three 
degrees of near wake. Assume the induced velocity is sought at the blade station, x = 9.5, and place 
the blade at azimuth \\i = 135 degrees. The accumulation of the shed wake’s influence is 

area under — — versus 9 
d9 

dv x 

The curve of — — versus wake age is illustrated with Fig. 29. This figure suggests an impending 
d9 

singularity as computations are made very close to the blade. In fact, the velocity induced at the 
blade station x = 9.5 becomes so large that a semi-log scale for the ordinate on Fig. 29 is helpful in 
capturing how rapidly the shed wake’s influence drops off as the blade moves away from the 
deposited shed wake. 

The common solution to the situation illustrated by Fig. 29 is to add a vortex core diameter, 
VCD, non-dimensionalized by rotor radius, R, to the L dimension of the shed-wake problem [refer 
back to the discussion surrounding Eq. (9.62)]. This solution is effective, as Fig. 39 shows. Of 
course, it is the integrated value 


/»vi/= 135 i 

1 180 av 


x=0.5 


180 


-132 


dO 


-dO 


that is more important, and this integrated result is shown versus vortex core diameter ratioed to 
radius in Fig. 31 for x = 0.5 and \|/ = 135 degrees. And even more important is the integrated value 
of induced velocity considering the wake extending all the way back to, say, 9 = -20 tt. The induced 
velocity at x = 0.5 (and with \|/ = 135 degrees) including this far wake is also shown on Fig. 3 1 . 


To examine the influence of the shed wake further, select, somewhat arbitrarily, the ratio of 
vortex core diameter to rotor radius equal to 0.015. Now the distance, L, between the vortex segment 
and rotor-blade station where induced velocity is sought, x, is rewritten as 


(0.102) L 3 =R 


[xv 2 + x 2 -2x(xv)cos(i|/-9)] 
+2p ( xv cos 9 - x cos v|/) (v|/ - 9) 

VCD V 


+ (v|/-9) (|A + kf pp ) + 


R 
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Mathcad’s built-in numerical integrator had absolutely no problem calculating the induced velocity 
at all radius stations and any azimuth. 


(0.103) 


v 


X,\| / 


"» 1 

C V 

( T poq fl ^ i 

I (2xv- x c -1)' 

x sin (v|/ - 0) + p (\|/ - 9) sin 0 


—00 

V A 1 COSU J 4 

1 1 r, 

(1-v) 2 

L L 3 


dOdxv 


The calculation could be performed even at x = xv because the distance, L, can never be smaller than 
VCD/R. Furthermore, because the shed wake’s circulation goes smoothly to zero at both the blade 
root and tip, there is no need to perform a fixed-wing coordinate transformation as was helpful for 
the trailed wake integration of both fixed and rotary wings. 


The next objective is to calculate the induced velocity for the sample problem’s geometry 
and operating condition. To begin the numerical integration, additional input is required, however. 
The additional input that satisfies rotor lift = 2,712 pounds at an advance ratio of 0.5 is that T 0 = 334 
and T i = -253.9. For convenience, the complete input for the shed-wake numerical integration is: 


R = 22 feet 

x c = 1/6 non-dimensional 
p = 0.002378 slugs per cubic feet 
VCD/R = 0.015 non-dimensional 
V t = 603.605 feet per second 
jlx = 0.5 non-dimensional 
Atpp = -0.03 non-dimensional 
r 0 = 334 square feet per second 
T i = -253.9 square feet per second 


With this infonnation as input, put Mathcad and its built-in integrator to work calculating Eq. 
(0.103) as follows: 
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INPUT 

p = 0.5 T o =334 r, =-253.9 R=22 x e =l/6 A. tpp =-0.03 VCD/R = 0.015 
Dimension integration. 

Select azimuth station, ¥ = variable for this problem 
Extent of wake age, WA = -207t 
Maximum number of radial stations, M = 90 

Range of radial stations where induced velocity is calculated, s = 0, 1 M - 1 


Then proceed with these calculations 


1 — x„ 


l + x„ 


K f 1 | 

a s = — I s + — =>To correspond to radius stations used in trailed- wake problem 


r fl + x c 

x. = — = 5 

s R l 2 


1 — x. 


cosa, =a-bcosa„ 


'V 


(F, cos9) 


L (2xv-x c -l) 


(!- x c) 


x sin ( \|/ - 0) + p( \|/ - 0) sin 0 

T 

I 5 

[xv 2 +x 2 -2x(xv)cos(\|/-0)j 
R-l +2p(xvcosO-xcosv|/)(\|/-0) 
fVCD ' 2 


d0dxv 


+ (\|/-0)"(p 2 +?t^ p ) 


l R 


x, =a-bcosa. 


The calculation of induced velocity at the lifting line — as created by the shed vortices from 
the lifting line — follows from the above scheme. Additionally, lift and horsepower distributions are 
calculated from 


( 0 . 104 ) 


al s , ¥ =[pV t (x s +gsin\|/)R]< (r 0 +r lS in ¥ ) i- (2x ; x ° ^ 

M-l 

L ¥ =X AL s v => Single-blade lift at input azimuth 

s=0 

AHP S , V = — v s. ¥ AL S jV => Running Horsepower 

M-l 

HP y = ^ AHP s => Single-blade Horsepower at input azimuth 

s=0 

where A , = ( 


) A s => Rutming Lift 
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The shed wake significantly contributes to the total induced velocity that the rotor blade sees, 
primarily over the mid-span portion of the blade. For example, the azimuthal variation in induced 
velocity at the three radius stations under examination (i.e., x = 0.25, 0.50, and 0.75) is provided by 
Figs. 32, 33, and 34. It is at the mid-span that the shed wake’s additional induced velocity is the 
greatest. 

The blade’s azimuthally varying lift for the balanced rotor is illustrated with Fig. 35. Note, in 
contrast to the unbalanced rotor shown on Fig. 24, the balanced rotor carries lift in the fore and aft 
quadrants of the revolution. This gives, in effect, a short wing span or low aspect ratio characteristic 
to the balanced rotor. As Fig. 36 shows, this concentration of balanced rotor lift in the fore and aft 
direction is accompanied by excessive induced horsepower when compared to the unbalanced rotor 
of Fig. 24. The average or steady induced horsepower for the balanced rotor is HP; = 76.4 hp versus 
50. 1 hp for the unbalanced rotor and versus 6. 13 hp for the ideal fixed wing. 

The requirement for a rotor to have zero rolling moment is clearly adverse to the induced 
power required to produce lift. The ratio of balanced rotor-induced horsepower to ideal fixed wing- 
induced horsepower is 76.4/6.13 = 12.4. But keep in mind that while an elliptical bound circulation 
is ideal for a fixed wing, it is not obviously true for a rotary wing. 
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V. Closing Remarks 


This elementary introduction to induced velocity has been presented assuming an elliptical 
bound circulation distributed along a lifting line to represent both the wing and the rotor. The 
assumption has been that 


T oc ^A n sinnp 

1 

and that Ai = 1 while A? through A* = 0. An extension of the fundamental equations provided by 
this note could easily be made using a full set of the Fourier series. This extension — if made — would 
completely generalize the lift distribution for the rotor just as is done in the case of fixed-wing 
theory. 


The lifting rotor’s spiraling vortex wake structure leads to very high induced power when 
compared to the ideal wing. For an advance ratio of 0.50, interference created by the spiraling rotor 
wake leads — just for a single blade — to induced power on the order of 10 times that of the wing 
when the comparison is made at wing span = rotor diameter and equal lift. While an elliptical bound 
circulation is known to be ideal for the fixed wing, it is quite probably not ideal for the rotary wing 
in high-speed forward flight. 

A single-bladed rotor and prescribed wake geometry have been selected for this rotary-wing 
introduction. This has been useful for an elementary discussion. However, the practical problem 
includes any number of blades and vortex wake structure that is free to deform based on 
fundamental principles. Furthermore, representing the rotor blade by a lifting line is quite 
unsatisfactory when the rotor lift distribution varies with time. Fortunately, advanced analyses 
coupled with powerful digital computers have given today’s rotorcraft engineers insight and 
practical answers to the effect of a rotor system’s wake upon the lifting surfaces that created the 
wake. 
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Figure 12 
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Figure 17 
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Figure 21 
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Figure 22 
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Figure 28 
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Figure 34 
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Figure 36 
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